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ABSTRACT 

Based on a carefully constructed sample of dwarf stars, a new optical-near infrared photo¬ 
metric calibration to estimate the metallicity of late-type K and early-to-mid-type M dwarfs 
is presented. The calibration sample has two parts; the first part includes 18 M dwarfs with 
metallicities determined by high-resolution spectroscopy and the second part contains 49 dwarfs 
with metallicities obtained through moderate-resolution spectra. By applying this calibration to 
a large sample of around 1.3 million M dwarfs from the Sloan Digital Sky Survey and the Two- 
Micron All Sky Survey, the metallicity distribution of this sample is determined and compared 
with those of previous studies. Using photometric parallaxes, the Galactic heights of M dwarfs 
in the large sample are also estimated. Our results show that stars farther from the Galactic 
plane, on average, have lower metallicity, which can be attributed to the age-metallicity relation. 
A scarcity of metal-poor dwarf stars in the metallicity distribution relative to the Simple Closed 
Box Model indicates the existence of the “M dwarf problem,” similar to the previously known 
G and K dwarf problems. Several more complicated Galactic chemical evolution models which 
have been proposed to resolve the G and K dwarf problems are tested and it is shown that these 
models could, to some extent, mitigate the M dwarf problem as well. 

Subject headings: galaxy: evolution - stars: late-type - stars: abundances - stars: fundamental 
parameters - techniques: photometric 


1. INTRODUCTION 

M dwarfs are the most numerous stars in the Galaxy, contributing about 70% of all stars by number 
(Reid & Gizis 1997, hereafter RG97). Their main-sequence (MS) lifetimes are much longer than the current 
age of the Universe and they can therefore be used as excellent tracers of Galactic structure and population 
as well as Galactic chemical, kinematical and dynamical evolution. Since the advent of deep, advanced 
surveys such as the Sloan Digital Sky Survey (SDSS, York et al. 2000) and the Two-Micron All Sky Survey 
(2MASS, Skrutskie et al. 2006), M dwarfs have been investigated in photometric and spectroscopic samples 
of unprecedented size, revolutionizing this area of astronomy. 

Clearly, a complete understanding of Galactic astronomy requires accurate knowledge of fundamental 
properties such as mass, radius, metallicity and temperature of these dwarfs. These, especially metallicity, 
however, have proven challenging to calibrate. Although accurate values of metallicity for M dwarfs can 
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directly be obtained by analyzing high-resolution spectra (e.g., Woolf & Wallerstein 2005, hereafter WW05), 
the development of alternative methods is needed. Since M dwarfs are among the intrinsically faintest stars , 
only a limited number of these stars are close enough for high-resolution spectroscopy (Woolf and Wallerstein 
2006). For this reason, efforts have been made to estimate the metallicity of M dwarfs based on spectral 
band indices through low-to-moderate-resolution spectra (Lepine et al. 2007, hereafter L07; Rojas-Ayala et 
al. 2010 & 2012; Terrien et al. 2012; Mann et al. 2013a, hereafter M13a; Mann et al. 2013b, hereafter M13b; 
Newton et al. 2014, hereafter N14; Mann et al. 2014, hereafter M14). 

Despite the significant progress in deriving empirical metallicity calibrations using moderate-resolution 
spectra within the last few years, there is still a need for simpler approaches to determine the metallicity of 
large numbers of M dwarfs. There have been Several attempts to do this using photometric properties. By 
employing a calibration sample of M dwarfs in M-I-FGK CPMSs, Bonfils et al. (2005, hereafter B05) derived 
an M-dwarf metallicity relation in terms of the A-band absolute magnitude (M^f) and V — K color with a 
dispersion of 0.2 dex. Johnson and Apps (2009, hereafter JA09) showed that while the relation of B05 could 
reasonably reproduce the metallicity of metal-poor M dwarfs, it underestimated the metallicities of their 
high metallicity stars by an average of 0.32 dex. They established an empirical model in which the distance 
of an M dwarf from the mean MS along in the {V — K) - plane was an indicator of its metallicity. 
Schlaufman & Laughlin (2010, hereafter SLIO) suggested that the empirical photometric calibrations of B05 
and JA09 systematically underestimated or overestimated metallicity at the extremes of their ranges. They 
improved upon those calibrations and determined the metallicity of an M dwarf using its distance from the 
mean MS along the V — K color (rather than M/f) in the {V — K) - M^k plane. By applying a new sample 
of M dwarfs in M-I-FGK CPMSs, Neves et al. (2012, hereafter N12) showed that the calibration of SLIO 
had a lower dispersion than those of B05 and AJ09, and slightly modified this relation by readjusting its 
coefficients using their own calibration sample. 

In order to use the photometric calibrations described above, the distances of M dwarfs are required; this 
limits their applications only to samples of stars whose parallaxes are available. A color-color diagram could 
provide a more efficient technique for determining the metallicity of M dwarfs, even those with unknown 
parallaxes. Several empirically determined metallicity calibrations through color-color diagrams have been 
derived in the last few years. Using their large spectroscopic sample of M dwarfs. West et al. (2011, hereafter 
Wll) obtained a two-dimensional fit which related a metallicity parameter (C, L07) to the g — r and r — z 
colors, with typical uncertainties of 10 - 20%. Bochanski et al. (2013) pointed out that the relation of Wll 
is limited only to near-solar metallicity M dwarfs. They introduced a quantity, <^(g-r) which measures the 
difference in g — r between a subdwarf and its solar metallicity counterpart as a function of r — z color. 
Most recently. West et al. (2014) used more than 20,000 M dwarfs with metallicities determined by the 
optimized spectroscopic calibration of M13b (tested using a sample of wide, low-mass binaries for which 
both components have an SDSS spectrum) and derived relations between the metallicity and the SDSS griz 
colors of M dwarfs. By employing an M-dwarf calibration sample in M-I-FGK CPMSs, Johnson et al. (2012, 
hereafter J12) developed a calibration which correlated the metallicity of an M dwarf with its distance to 
the MS along the J — K color in the (V — K) - (J — K) plane, with an RMSE = 0.15 dex. N14 assembled a 
sample of 447 M dwarfs with metallicities calculated by their own spectroscopic calibration and found that 
the J — K color of an M dwarf is the best single-color diagnostic of its metallicity. They then derived a 
metallicity relation as a function of the distance of an M dwarf to the MS along J — K color in the {J — K) 
- {H — K) color-color diagram, with a multiple correlation coefficient (R^p) of 0.92. 

Our main goal in this study is to derive an optical-NIR photometric metallicity calibration which can 
readily be applied to large numbers of stars without the need for parallaxes, moderate-to-high resolution 
spectra or time-intensive calculations. This will allow us to study statistically the metallicity distribution of 
the local Galactic disk and test models of Galactic chemical evolution (GGE). The observations, sampling 
process and M-dwarf selection are briefly described in section 2. Our metallicity calibration sample and 
the best-fit relation for estimating M-dwarf metallicity are presented in Section 3. The Galactic height {Z 
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- height) distribution is investigated in Section 4. Section 5 is devoted to the application of the metallicity 
relation to our large sample and the comparison of the resulting metallicity distribution with those from 
other studies. The statistical relation between the metallicity and Z - height of M dwarfs is also discussed 
in this Section. In Section 6, the Simple Closed Box Model (SCBM, Schmidt 1963) as well as a few more 
realistic GCE models are compared with our M-dwarf metallicity distribution. 


2. OBSERVATIONS AND M DWARF SAMPLING 
2.1. Matched SDSS And 2MASS Data 

The SDSS is one of the most extensive surveys in astronomy. During its operations, it obtained multi¬ 
color images which covered more than a quarter of the sky. Multi-band photometry was collected using 
a 2.5-m wide-angle optical telescope at Apache Point Observatory, in New Mexico. The camera included 
thirty CCD chips each with 2048x2048 pixels (with a total of 120 Megapixels). These chips were organized 
in five columns of six chips per column. Each column had a different optical filter bandpass, designated 

u, g, r, i, z (Fukugita et al. 1996), with average wavelengths of 355.1, 468.6, 616.5, 748.1 and 893.1 nm 
and with 95% completeness of point sources in typical seeing to magnitudes of 22.0, 22.2, 22.2, 21.3, and 
20.5, respectively (Gunn et al. 1998). The SDSS imaging covers several thousand square degrees of sky, 
and over this region, photometric calibrations achieved an accuracy of « 0.01 - 0.03 mag in ugriz (Ivezic 

al. 2004, 2007; Adelman-McCarthy et al. 2006; Tucker et al. 2006; Padmanabhan et al. 2008). 

Our calibration uses both optical (the SDSS g) and NIR (the 2MASS JHK) broadband filter and we 
cross-matched the SDSS and 2MASS data for collecting our sample. The 2MASS was a survey of the whole 
sky in the three IR broad-band filters J, H and K, with average wavelengths of 1.25 /xm, 1.65 ^m and 
2.17 /xm, respectively. The measurements were made from 1997 to 2001 using two highly-automated 1.3-m 
telescopes, one at Mt. Hopkins, Arizona for the Northern Hemisphere observations and one at Gerro Tololo 
Inter-American Observatory, Chile for the southern hemisphere data. Each telescope had a three-channel 
camera, each containing a 256 x 256 array of infrared detectors to scan the sky simultaneously in the three 
filters. Point sources brighter than about 1 mJy, with S/N greater than 10 in each band were characterized 
and compiled in a separate catalog, the 2MASS Point Source Catalog (2MASS-PSC, Cutri et al. 2003). 
These bright sources (< 13 mag) generally have 1-sigma photometric uncertainty of < 0.03 mag (Skrutskie 
et al. 2006). 

We used the SQL Search tool in the SkyServer DR9 to collect an M-dwarf sample suitable for our study. 
The stars were taken from the SDSS DR9 photometric catalog (identified by PhotoObj, Ahn et al. 2012) 
with the object class specihed by Type = 6 (for stars). Each object in the SDSS DR9 catalog has a unique 
SDSS identifier, called ObjID, which is also included in the 2MASS-PSC for objects in the regions that 
overlapped with the SDSS. We therefore selected those stars with the same ObjIDs in the two catalogs. 


2.2. Extinction Correction 

To minimize extinction, we chose stars with high Galactic latitude b > 50°. It has long been shown 
that, on average, extinction decreases at higher latitudes, due to the low column densities found along these 
lines of sight (Zagury 2006;Larson and Whittet 2005). 

We employed the dust maps of Schlegel et al. (1998, hereafter S98) to correct the stellar photometry 
for Galactic extinction. Gurrently, these maps provide the most comprehensive dust data of the Milky Way 
on a large scale, including two two-dimensional full-sky maps (one for the northern and one for the southern 
Galactic hemisphere) with total line-of-sight dust column densities determined from far-infrared (100 and 
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240 /im) emission data. These maps have a resolution of 6.1 arcmin and are shown to predict reddening 
within 16% (S98). We also used the relative extinctions given in Table 6 of S98 to convert the extinction 
corrections in the Wband magnitude to SDSS-2MASS filter bands. While the maps mentioned above provide 
an efficient tool for estimating dust extinction, they refer to the total Galactic extinction along the line of 
sight and so may overestimate the true extinction to nearby stars (Covey et al. 2007 hereafter C07, Jones 
et al. 2011 hereafter Jll). Most of objects in the Galaxy lie behind only a portion of the Milky Way’s total 
dust column and thus are likely attenuated and reddened by only a fraction of the total dust column. 

Using spectra of more than 56,000 M dwarfs from the SDSS, Jll created a high-latitude, three- 
dimensional extinction map of the local Galaxy. In their technique, spectra from stars in the SDSS DR7 
dwarf sample along low-extinction lines-of-sight were compared with other SDSS M-dwarf spectra for deriv¬ 
ing distances and accurate lines-of-sight extinction. The three-dimensional map is most appropriate for stars 
within around 500 pc from the Galactic plane (Z < 500 pc). As mentioned in Jll, it can safely be assumed 
that stars with Z - heights greater than 500 pc are essentially behind the entire dust column and the maps 
of S98 are more useful for such stars. 


2.3. Clean Photometry 

Due to partial overlaps between neighboring SDSS images, an object might be observed more than 
once. The best observation of an object is called PRIMARY detection and the other observations, if any, 
are assigned as SEGONDARY detections. To avoid duplication, we chose only PRIMARY observations by 
setting the variable “mode” to 1. We also used point-spread function (PSF) magnitudes for SDSS data since 
these optimally measure the total fluxes for stars in this study. 

By equating the flag CLEAN to 1, we selected those stars whose SDSS magnitudes has passed appro¬ 
priate standards of clean photometry (the flags: NOPROFILE, PEAKCENTER, NOTCHECKED, PSF- 
FLUX-INTERP, SATURATED, BAD-COUNTS-ERROR, DEBLEND-NOPEAK and INTERP-CENTER 
are not set). We also set the 2MASS read flag to “222”, blend flag to “111”, and contamination/confusion 
flag to “000” to select only those stars which have unsaturated, unblended and uncontaminated photometry 
in the three IR bands {J, H , and K). 


2.4. Dwarf-Giant Separation 

We selected stars which fall within typical color ranges for late-type K and M dwarfs, r — i > 0.47 
and i — z > 0.25 (Wll). There possibly may be contamination by giants in any color-selected sample of M 
dwarfs, however, which needs to be addressed in statistical studies of the Galaxy. Bessell & Brett (1988, 
hereafter BB88) found in the JHK, there is a clear bifurcation between M giants and M dwarfs, beginning 
where TiO bands appear or the M-dwarf sequence starts. By applying the NIR color ranges typical of M 
giants (BB88; C07), we found that around 1% of stars in our SDSS-2MASS sample are giants, consistent 
with other studies such as Covey et al. (2008, hereafter C08) with less than 2% and Wll with around 0.5% 
giant contamination. 

It should be mentioned that this approach for separating giants from dwarfs is not accurate for K-type 
stars (BB88) as they overlap in JHK color space. A better way to remove giants employs reduced proper 
motions of stars as described in M13b. By using a sample of stars with available proper motions, giants 
could be separated from dwarfs with higher accuracy. 
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2.5. Photometric Parallax 

The distance and Z - height of our stars were estimated by means of a photometric parallax technique 
developed by Bochanski et al. (2010, hereafter BIO). More specifically, the absolute magnitude in the r 
band, M,., of low-mass stars can be obtained from a (J — K)— relation (the first equation of Table 4 in 

BIO) which is valid for the color range 0.50 < {r — z) < 4.53. We therefore selected only those stars which 
fall in this color range. 

After meeting all the criteria above, as well as the color cuts applicable to our metallicity calibration 
outlined in Section 3, a sample of 1,330,179 M dwarfs was selected. 


3. METALLICITY CALIBRATION 

To calibrate metallicity using photometry, we collected a calibration sample of M dwarfs with reliable 
metallicity and photometry. Since the dust maps of S98 were used to remove Galactic extinction for our large 
sample, the metallicity calibrations to be used for this sample must be based on calibrators with photometry 
corrected in the same way. For this reason, we corrected the photometry of all stars in our calibration sample 
using these map^ 

There are nearby M dwarfs which are members of binary systems, having FGK-dwarf primaries with 
metallicities determined through high-resolution spectra. We identified 14 such M dwarfs with unsaturated 
SDSS g and unsaturated 2MASS JHK magnitudes while for most of these stars the photometry in the r, i or 
z bands are saturated. We used the SDSS_Phot Flag to check the saturation status of stars; the magnitude 
in a SDSS filter band is unsaturated if the SATURATED flag for this band is not set. To be certain of 
unsaturated JHK photometry, we chose only those stars whose 2MASS read flag is “222”. 

In order to add more stars with reliable metallicities to our calibration sample, we refer to the work of 
Dahab & Straus^ (personal communication) who showed that by filtering the saturated objects in SDSS 
that have certain flags set, the photometry of remaining saturated objects are usable. The stellar locus 
on a color-color diagram can be identified by the locations at which photometrically clean stars appear. 
By comparing color-color diagrams of stars with unsaturated and saturated photometry, Dahab & Strauss 
demonstrated that the saturated PSF magnitudes in the ugri bands for which none of the flags EDGE, 
INTERP-CENTER and PSF-FLUX-INTERP are set leave the color-color diagrams with almost no oulier. 
They also pointed out that the g magnitudes of these saturated stars have a minimum of about 13.4 mag 
and the number of stars fainter than this minimum increases exponentially. This indicates that saturated 
stars with g < 13.4 mag are problematic and must be excluded from photometric studies. 

An object which is too close to the edge of an image is flagged EDGE in the SDSS photometry. Among 
PRIMARY objects, only large extended objects should be flagged EDGE. Therefore, for point sources, there 
is no need for concern about the PRIMARY objects having the EDGE flag set. As a result, rather reliable 
photometry is possible after filtering out the saturated objects in the SDSS whose INTERP-GENTER and 
PSF-FLUX-INTERP flags are set. In this way, we found three M dwarfs in M-I-FGK CPMSs, which have 
unsaturated photometry in the JHK bands but whose g magnitudes are saturated while none of the two flags 
are set. Moreover, the g magnitudes of these three stars are fainter (by around 4 mag) than the lower limit 
mentioned above, and we thus call them “weakly saturated” magnitudes. According to Dahab & Strauss, 


dust map may make an offset in the extinction value of stars, specifically for each filter band. To minimize this effect 
on metallicity, the same offset must be included in the extinction of stars in both the calibration sample and large sample. For 
this reason, the same dust maps were applied for both samples. 

^Unpublished undergraduate Honor’s thesis of W.E. Dahab under supervision of M.A. Strauss, Princeton University. 
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Fig. 1.— The {g — K) - (J — K) color-color diagram for the 18 M dwarfs with metallicities determined 
by high-resolution spectroscopy. The metallicity values are color-coded: stars with [Fe/H] > -1-0.15 dex are 
plotted in red, with —0.1 < [Fe/H] < -1-0.15 dex in green, with —0.4 < [Fe/H] < —0.1 dex in purple and with 
[Fe/H] < —0.4 dex in blue. The spectral types are symbol-coded: spectral types between around MO and 
M2 are shown by closed triangles, around M3 by closed diamonds, between around M4 and M5 by closed 
squares, and around M6 by closed circles. 


we can be certain that these magnitudes are reliable enough for our study. 

The sample of M dwarfs in CPMPs now includes 17 (14-1-3) stars having metallicities determined through 
high-resolution spectra. Moreover, we found the measeured photometry and metallicity of LHS 3084 reliable 
and added this to the sample, increasing the number of stars to 18. Although LHS 3084 is not in a CPMP, its 
metallicity was determined by direct measurements on its high-resolution spectrum, not that of a primary. 
The astrometry, extinction-corrected photometry and saturation status in the g band for the 18 M dwarfs 
are shown in Table 1. The spectral type of these M dwarfs as well as the spectral type and metallicity of 
the primaries are given in Table 2. This Table also includes the references from which the metallicities and 
spectral types are taken. 

The 18 M dwarfs described above form the foundation of our metallicity calibration. We investigated the 
relation between the metallicities and all possible colors involving g, J, H and K magnitudes for these stars. 
It was found that the J—K color is the best metallicity indicator which is in agreement with the previous 
studies of J12 and N14. There are deep potassium (K) and iron hydride (FeH) absorption features as well 
as dozens of shallower metal lines in the J band spectra of M dwarfs. Consequently (as J12 suggested), M 
dwarfs with higher metallicity have preferentially suppressed J band fluxes as compared to A-band spectra 
where there are only a few relatively shallow Na and Ca lines as the prominent absorption features. This 
causes the J — K color of the metal-rich M dwarfs to be redder. 

In addition, we found that among all possible color-color diagrams, the {g — K) - {J — K) diagram can 
reliably separate metal-poor from metal-rich M dwarfs. Figure 1 shows such a diagram for the 18 M dwarfs. 
The metallicity values are coded by colors and the spectral types are coded by symbols provided in the 
caption. Overall, for a given spectral type range, metal-poor stars have bluer J — K than metal-rich ones, 
and clearly the J — K color is a better diagnostic for metallicity than g — K. It can also be seen that the 
g — K color is a better indicator for spectral type than J — K^ and for a given metallicity range, earlier-type 
M dwarfs have bluer g — K than stars of later spectral type. The location of a star in a color-color diagram 
theoretically depends on its fundamental properties such as metallicity, spectral type (or temperature) and 
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Fig. 2.— The {g — K) - (J — K) color-color diagram for the 67 dwarfs in the calibration sample. The 
metallicity values are color-coded as described in Figure 2. The spectral types are symbol-coded: spectral 
types between around K6 and K7 are shown by hexagrams, between around MO and M2 are shown by 
triangles, around M3 by diamonds, between around M4 and M5 by squares, and around M6 by circles. The 
18 stars from Figure 1 are depicted by closed symbols. 


surface gravitj0 However, as an approximation, we assumed log((7) is constant for all M dwarfs (« 5) and 
has the same effect on the colors of these dwarf stars. Generally, different colors have different sensitivities 
to metallicity and spectral type. In our case, while the J — K could be a measure of metallicity, the g — K 
is not sensitive enough to metallicity to separate metal-rich from metal-poor M dwarfs adequately. On the 
other hand, the g — K color can distinguish different spectral types more effectively than the J — K color. 

Although the 18 M dwarfs provided us with a sample of robust calibrators, these stars by themselves are 
insufficient for covering the entire {g — K) - (J—K) plane, however, as evidenced by gaps, particularly for [g — 
K) values between 7 and 8. If we could add more dwarfs with trustworthy photometry and metallicity, there 
would be the same trends as those of the 18 M dwarfs: for a given spectral type range, more metal-poor stars 
should have smaller values of J — AT, and for a given metallicity range, earlier-type stars should have smaller 
values of g — aQ This defines a rule for selecting the rest of stars in our calibration sample, we collected 
55 early-to-mid-type M dwarfs and late-type K dwarfs with spectroscopically determined metallicities from 
moderate-resolution spectra. Their selection initially was by eye, based on their locations on the {g — K) 
- {J — K) diagram with respect to the original 18 M dwarfs. We chose those dwarfs which approximately 
followed the rule described above and this enabled us to exclude some outliers from the sample in the first 
step of the selection. Assuming accurate photometry, large uncertainties in metallicity and spectral type 
could show a star in a metallicity-spectral type category different from its real one, causing it to significantly 
deviate from the rule. 

We then obtained the best fit of a low-order polynomial between metallicity and g — K and J — K colors 
for the whole sample (18-1-55=73 stars) as follows: 


®These are the three main properties which play important roles in constructing model atmospheres and synthetic spectra 
of stars. 

^Similar to this can, to some extent, be perceived from the work of N14 (see Figure 21 of their paper) on the (J — K) - 
(H — K) diagram. 




\FelH\ = Cl + C2(g - X) + C^{J -K)f - Kf + C^{J - Kf + - K){J - K) (1) 

where Q are constant coefficients. This relation is of the same order (second order) as that of J12 who used 
an optical-NIR calibration with a filter set similar to ours. The only difference is that we employed the g 
magnitude instead of the V magnitude which was used in the calibration of J12. We rejected 6 outliers for 
lying more than 2.5-sigma from the relation, leaving a final calibration sample of 67 stars. The astrometry, 
photometry, saturation status in the g band, spectral type, metallicity and related references for the rest of 
our calibrators are listed in Tables 3 and 4. The {g — K) - {J — K) color-color diagram of the whole sample 
is presented in Figure 2. 

The coefficients of the final best fit for the 68 dwarfs in the calibration sample are: 

Ci = {-14.2959,0.0519,29.5926,-0.0529,-17.6762,0.7032} 

with an RMSE = 0.077 dex and 7?^^ = 0.90 (which are an improvement upon the values for the sample of 
the original 18 M dwarfs: RMSE=0.090 dex and R'^p = 0.086), yielding elliptical isometallicity contours. 
It should be remarked that this metallicity calibration is applicable for stars of spectral types between K6 
and M6.5 with -0.73 < [Fe/H] < -bO.3 dex, 3.37 < g - K < 8.46 and 0.71 < J - K < 1.01. Figure 3 
shows a comparison between the values of metallicity obtained by Equation (1) with coefficients above and 
those taken from other studies for the same 67 calibrators. We should mention that to have a more accurate 
metallicity calibration, a careful investigation of surface gravity is also required. 

In the following sections, this metallicity calibration is applied to our large sample and the resulting 
metallicity distribution is studied in more detail. 



- 0.4 - 0.2 0 

Metallicity - This Work 


Fig. 3.— The metallicities derived by other studies versus metallicities calculated by Equation (1) 


4. Z - HEIGHT DISTRIBUTION 

To study the relation between metallicity and Z - height, it is necessary to estimate the distance of 
stars in our sample. We applied Equation (1) to the sample of 1,330,179 M dwarfs described in Section 2, 
accepting only those stars with the metallicity range within which this equation is valid. This left a sample 
of 1,298,972 stars: only ~ 1% of stars have metallicities outside this range; the rest satisfy all requirements 
outlined in Section 2. We also limited our study to stars with Z < 2000 pc, leading to a final sample of 
1,298,721 M dwarfs (hereafter, “SampleMetal”, or briefly “SM”), having a median extinction in the r band 
of Ar = 0.058 mag. 
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We estimated the Z - height of stars in SM using the photometric parallax described above (Subsection 
2.5). It should be remarked that these Z - heights are the vertical distances from the plane passing through 
the Sun. Since all the selected M dwarfs lie above this plane {b > 50°), we added the Z - height of the Sun, 
Z:^ = 20 pc (Juric et al. 2008), to those derived from the photometric parallax to obtain the distances from 
the central plane of the disk. 

The Z - height distribution of the sample (Figure 4) shows a peak at Z ~ 350 pc, indicating that around 
94% of the stars are within the Galactic thin disk, i.e., Z < 1000 pc (Gilmore & Reid 1983). 



Fig. 4.— The Z - height distribution of stars in SM from the observed data 



Fig. 5.— The metallicity distribution of stars in SM from the observed data (dashed histogram) and volume- 
corrected data (solid histogram) 


5. METALLICITY DISTRIBUTION 

The metallicity distribution of stars from SM having —0.7 < [Fe/H] < -1-0.3 is shown as a dashed 
histogram in Figure 5. As mentioned above, only those stars with b > 50° were selected. This means 
that stars outside an imaginary cone, perpendicular to the Galactic plane with an opening angle of 40° 
and apex at Earth, were excluded. As a result, stars with higher vertical distances are distributed through 
larger volumes of space. In order to remove the effect of this bias on metallicity distribution, it is necessary 
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to perform a volume correction. To this end, the stars in SM were divided into 20 bins of 100 pc in Z - 
height and the weighting factor proportional to the inverse of the volume of each bin was calculated. The 
metallicity distributions of these bins were obtained separately and multiplied by the corresponding weighting 
factors. All these weighted distributions were then added together as the total volume-corrected distribution, 
illustrated in Figure 5 (the solid histogram). This distribution shows a maximum at around -1-0.1 dex, and 
a mean metallicity [Fe/H] « —0.04 dex. Note that this mean metallicity is not expected to represent the 
true mean value of the local Galactic disk since those stars with metallicity outside the metallicity range of 
our calibration were removed from the sample. In order to obtain a more meaningful mean metallicity, a 
calibration valid over a wider range of metallicity is necessary. 


5.1. Metallicity and Z - height Relation 

The relation between metallicity and distance Z above the Galactic plane can be investigated by deter¬ 
mining the metallicity distribution of stars as a function oi Z - height. We divided a subsample of our stars 
with Z < 1200 pc into 12 bins of equal AZ = 100 pc. The metallicity distribution associated with each bin 
is shown in Figures 6 and 7. There is a clear shift towards lower metallicities as Z - height increases.; the 
fraction of metal-rich stars slowly decrease and the distributions become more metal-poor with increasing 
height. The mean metallicity decreases (by around 0.19 dex) from 0 to 1200 pc which is consistent with the 
age-metallicity-height relation. Studies of Galactic evolution have demonstrated that stars which formed at 
earlier times in the Galaxy’s history generally have lower metallicities in general and are, on average, farther 
from the Galactic plane. Stellar systems are formed from interstellar gas and, at the end of their lifetimes, 
they may return a substantial fraction of their initial masses enriched with metals to the interstellar medium. 
As a result, succeeding generations of stars become more metal-rich than their predecessors, leading to an 
age-metallicity relation. The pioneering work on this was done by Twarog (1980a, 1980b) who derived the 
age-metallicity relation for the disk in the neighborhood of the sun using a large sample of southern F dwarfs. 
They found that the mean metallicity of the Galactic disk increased by about a factor of five between 12 
and 5 billion years ago and has changed only slightly since then. A recent evidence of this kind was offered 
by Gasagrande et al. (2011, hereafter Cll) who determined the metallicity distributions of three samples of 
solar-type stars with different ranges of age (Figure 16 of their paper). They showed that young stars (ages 
< 1 Gyr) have a quite narrow distribution around higher values of metallicity, whereas intermediate-age 
(between 1 and 5 Gyr) and old-age (> 5 Gyr) stars present broader distributions, significantly extending to 
lower metallicities. 

In addition to the age-metallicity relation, there is also a relation between age and velocity dispersion. 
Observations have demonstrated that old stellar populations in the disk have larger velocity dispersions. 
Based on the sample of Twarog above, Garlberg et al. (1985) examined the age-velocity dispersion relation 
for stars in the solar neighborhood. They found an increase in velocity dispersion to a maximum at about 6 
Gyr, thereafter, the trend remained roughly constanlj^ The term “dynamical heating” is commonly applied 
to all processes that cause an increase in velocity dispersion with age. Different mechanisms are responsible 
for vertical dynamical heating (e.g., Nordstrom 2008). Whatever the mechanism(s), stars dissipate from the 
Galactic plane, increasing their vertical distance from the plane in the course of time. One therefore expects 
that older, more dynamically heated stars should be found, on average, at larger Z - heights than the younger 
stars. As a consequence, stars farther from the Galactic plane should be, on average, more metal-poor than 
those closer to the plane. 


®This behavior is consistent with a simple model for the local-disk formation and evolution in which growth occurs at a 
uniform rate with no initial disk, and the random velocities of stars are increased by a mechanism which depends on the total 
surface density of the disk, such as spiral waves. 



Fraction Fraction Fraction 
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Metallicity 



Fig. 6.— The metallicity distributions of M dwarfs in the large sample with different Z - height ranges 
between 0 and 600 pc. The Z - height range and number of stars for each distribution are labeled on the 
corresponding plot. 
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Metallicity 



Metallicity 


Fig. 7.— The metallicity distributions of M dwarfs in the large sample with different Z - height ranges 
between 600 and 1200 pc. The Z - height range and number of stars for each distribution are labeled on the 
corresponding plot. 
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Using a metal-sensitive ratio, (CaH2 -|- CaH3)/Ti05, from Lepine et al. 2003 as a proxy for M-dwarf 
metallicity, West et al. (2008) demonstrated a decrease in metallicity as a function of Z - height (up to 1000 
pc) which implies that stars more distant from the Galactic plane have lower metallicities. In a simulation 
by West et al. (2006; 2008), a simple one-dimensional model of the thin disk was applied to investigate the 
vertical motions and positions of M dwarfs over the Galaxy’s lifetime. The model also showed a decrease 
in activity fractions (as traced by Ha emission) as a function of Z - height for all M-dwarf spectral types, 
suggesting a decline in the magnetic activity of older M dwarfs. All results from this ID dynamical simulation 
are consistent with observed activity and velocity trends. 


5.2. Comparison With Previous Studies 

Figure 8 indicates a comparison between the metallicity distribution (volume-corrected) from this work 
(solid line) and that from WW12 who fobtained the metallicity distribution of 4141 M dwarfs in the spectro¬ 
scopic SDSS catalog (dashed line). As shown in the Figure, there is a slight offset between the peaks of the 
two M-dwarf distributions. It appears, however, that our sample includes a larger fraction of metal-poor stars 
than that of WW12. Moreover, the distribution of WW12 has a steeper slope towards lower metallicities 
than ours. 



Fig. 8.— The metallicity distribution of the volume-corrected sample from this work (solid histogram) and 
the distribution from WW12 (dashed histogram) 

Possible reasons for the discrepancy between these two distributions are as follows. First, the distri¬ 
butions were obtained by different methods. To estimate the metallicity of their sample, WW12 derived a 
linear relation (Equation (1) in their paper) between [Fe/H] and the metallicity-sensitive parameter C defined 
in L07. This relation was applicable to M dwarfs with 3500 < Tes < 4000 K and —1.5 < [Fe/H] < -1-0.05 
dex. Dhital et al. (2012) recalibrated the definition of ( and showed that this new ( could be a significantly 
better indicator of metallicity for early M-type dwarfs (between MO and M3). M13a tested the parameter 
defined in L07 using their own sample and found that it is not only sensitive to metallicity (and temper¬ 
ature) but also to some other stellar characteristics such as activity and surface gravity. This could lead to 
an incorrect identification of some metal-poor stars as near-solar metallicity stars. They pointed out that 
although C correlates well with [Fe/H] for supersolar metallicities, it does not always diagnose metal-poor M 
dwarfs correctly and may identify such stars as more metal-rich than they are. Further, Lepine et al. (2013) 
remarked that the parameter C defined in L07 overestimated the metallicities of early-type M dwarfs while 
the parameter defined in Dhital et al. (2012) underestimated the metallicities of these stars, leading them 
to the redefinition of the index (. Due to the overestimation of the metallicity of metal-poor M dwarfs by 
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the parameter ^ from L07, we could expect a bias toward metal-rich stars in the metallicity distribution of 
WW12, as seen in Figure 8. 





Fig. 9.— The metallicity distribution of the volume-corrected sample of M dwarfs from this work (solid black 
histogram) and those of G dwarfs (dashed green histogram) and K dwarfs (dashed red histogram) from S12 
for different height ranges. 

Second, there is an important difference between our sample and that of WW12. As shown in Section 
4, our sample is taken mostly from the Galactic thin disk, and a large fraction of its stars are in the solar 
neighborhood with Z < 500 pc. On the other hand, a significant portion of stars in the sample of WW12 lies 
outside the thin disk and is not representative of the local neighborhood: the WW12’s sample is expected 
to have more low-metallicity stars than ours, but that is not what is seen in Figure 7. 

It should be noted that the distribution of WW12 in Figure 8 is volume-corrected due to the metallicity- 
dependent volume coverage. It should be mentioned that this is different from the volume correction (owing 
to only volume effects) which we made above. WW12 argued that since low-metallicity main sequence 
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stars (subdwarfs) are less luminous than higher-metallicity stars of the same temperature or spectral type, 
higher metallicity stars should be overrepresented in magnitude-limited samples. In other words, metal- 
poor M dwarfs must be, on average, closer than the metal-rich ones in such samples. Using a sample of 
stars from Woolf et al. (2009) for which parallax data were available, WW12 determined the luminosity 
variation with metallicity for M dwarfs in the temperature range 3500 < Tgff < 4000 K. Based on this 
variation, they calculated volume-correction factors for different values of [Fe/H] (Table 1 in their paper) 
and implemented a volume correction for their metallicity distribution. Accordingly, to have a more reliable 
metallicity distribution, we need to determine the metallicity-dependent correction factors appropriate to our 
magnitude-limited sample. Nevertheless, we would not expect these corrections to change the location of the 
peak and the slope of the low-metallicity tail significantly (Figure 2 of WW12). Therefore, the uncorrected 
metallicity distribution can give us, to some extent, important information about the local Galactic disk. 

Figure 9 compares our M-dwarf metallicity distributions with those of G and K dwarfs in the SEGUE 
survey from Schlesinger et al. (2012, hereafter S12) for three different Galactic height ranges. As noted in 
S12, the G- and K- dwarf metallicity distributions have more metal-poor stars as height increases. It can be 
seen that the M-dwarf distribution exhibits the same trend, which implies a similar star formation history 
as G and K dwarfs. There is also a discrepancy between the distributions of different spectral types; the 
cooler spectral type, the more metal rich the distribution includes. 

JA09 argued that the M and FGK dwarfs in the local Galactic thin disk should have the same metal¬ 
licity distributions and that there should not be a systematic offset between M- and FGK-dwarf metallicity 
distributions. They related the offset between the metallicity distributions of M dwarfs and FGK dwarfs 
in some published papers to the reliability of the methods by which the metallicity of M dwarfs were de¬ 
termined. Accordingly, they adopted the mean M dwarf metallicity of their sample to be the same as the 
mean metallicity of a volume-limited sample of FGK dwarfs (~ —0.05 dex) from VF05. SLIO improved 
upon JA09 and remarked that for a fair comparison between samples of M dwarfs and FGK dwarfs in the 
local Galactic disk, the samples must not only be volume corrected but also have equivalent kinematics. By 
selecting two volume-completed, kinematically-matched samples of FG and M dwarfs, SLIO derived a mean 
metallicity of around —0.14 dex for M dwarfs in the solar neighborhood. It should be mentioned that we 
compared our distribution to those of previous works without any consideration of kinematics of the samples 
under study. To make proper comparisons between the metallicity distributions of dwarf stars with different 
spectral types, volume-corrected samples with equivalent kinematics have to be considered. 

On the other hand, the offset between FGK- and M-dwarf distributions might be explained in other 
ways. For example, similar to S12, one may expect that since M dwarfs have longer lifetimes than G and 
K dwarfs, there is a possibility that more metal-rich G and K dwarfs have evolved off the main sequence, 
causing a metallicity bias in G- and K-dwarf distributions against high-metallicity values which has not 
happened for cooler M dwarfs. Furthermore, if more lower-mass dwarfs are born than higher-mass stars as 
the metallicity of interstellar gas increases, there should be fewer metal-poor M dwarfs than G and K dwarfs. 

While there have been discrepancies between the metallicity distributions of different spectral types 
from various studies and some effort has been made to explain them, we can mention the work of GlI whose 
metallicity distribution of solar type stars (obtained by a photometric method) has a peak around [Fe/H] ~ 
0, close to those of the M-dwarf distributions shown in Figure 8. This is in agreement with the statement 
of JA09 who pointed out there should not be any offset of different metallicity distributions. To reach an 
accurate conclusion, more careful investigations are needed. 
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6. THE Galactic Chemical Evolution (GCE) MODELS 
6.1. The Simple Closed Box Model (SCBM) 

The simplest model of GCE, the SCBM, is based on three assumptions: 

1. The system is isolated and there is no mass inflow or outflow, 

2. There are two types of stars; those that live forever and those that die right after their birth, indicating 
instantaneous recycling, 

3. The system initially starts off with entirely metal-free gas and eventually ends up full of stars. 

Studies of the chemical evolution of local G dwarfs (e.g., van den Bergh 1962; Pagel & Patchett 1975; 
Wyse & Gilmore 1995) have shown the existence of the “G dwarf problem”, which indicates that the SCBM 
predicts many more low metallicity G dwarfs than are observed. Similarly, the “K-dwarf problem”, a 
paucity of metal-poor K dwarfs relative to the prediction of the SCBM, has also been observed (e.g., Casuso 
& Beckman 2004, hereafter CB04). WW12 were the first to recognize the “M-dwarf problem” (as shown in 
Figure 2 of their paper); the number of low metallicity M dwarfs is insufficient to match the SCBM, as G 
and K dwarfs. 

To test the SCBM, we compare our metallicity distribution with that predicted by the model. It can be 
shown (e.g., Pagel 2009 hereafter P09; Mo et al. 2010) that the mass in stars (Ms) with metallicities between 
Z and Z -|- dZ is given by 


dMs oc exp(— 



( 2 ) 


where p is the metal yield, defined as the mass of newly synthesized and ejected metals per unit mass 

z 

abundances between [Fe/H] and [Fe/H]-|-d[Fe/H] can be estimated as 


permanently locked up in the stars. Considering [Fe/H] ~ [M/H] ~ log( 2 ^), the mass of stars with metal 


dMs oc Zexp(— 


p)d[Fe/H], 


which can be rewritten as 


dMs 

cllog(Z) 


oc Zexp(-) 

P 


( 3 ) 


( 4 ) 


The proportionality (4) is the metallicity distribution function predicted by the SCBM. The histogram 
in Figure 10 depicts our volume-corrected metallicity distribution. The black curve in this Figure shows 
the normalized distribution from the SCBM with p=0.025 which has a maximum at the same metallicity 
([Fe/H]=-|-0.1 dex) as our histogram. 

Obviously, the SCBM cannot reproduce the low metallicity tail of the distribution, suggesting the M 
dwarf Problem. It is clear that one or more of the assumptions in the SCBM, such as instantaneous recycling, 
the system’s isolation without gas flow or the assumption of a zero initial metallicity are unrealistic. The G 
and K dwarf problems have motivated a generation of more complicated GCE models. We will examine a 
few such models and test them using our M-dwarf metallicity distribution in the following subsections. 
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Fig. 10.— The comparison between the metallicity distribution of M dwarfs from this study (the histogram) 
and those from the SCBM with different initial metallicities: Zq = 0 (the black curve), 0.002 (the red curve), 
and 0.003 (the blue curve). 


6.2. The SCBM With Pre-enrichment 

One of the solutions proposed for the G dwarf problem is to assume an initial metallicity of Zg 7 ^ 0 
(Samland & Hensler 1996). If, for some reason, there is a finite initial abundance of order [Fe/H]o « —1 dex, 
a good fit to the distribution is obtained at the expense of rejecting some stars with low metallicity. In this 
case, all metallicities (Z) in Equation (4) are replaced by Z — Zg yielding: 


dMs 

dlog(Z) 


a (Z - Zg)exp(--^^--^) 


(5) 


The normalized distributions calculated by the SCBM with two selected initial metal abundances of Zg 
= 0.002 (the red curve) and 0.003 (the blue curve) are shown in Figure 10. We set the yields for the red 
and blue curves to p = 0.023 and 0.022, respectively, to have a maximum at [Fe/H]=+0.1 (the same as our 
distribution). As can be seen, the distributions of the SCBM with a pre-enrichment give a better fit to our 
measured metallicity distribution than that without any initial abundances (the black curve). 

There have been several hypotheses to justify such a pre-enrichment. For example, this might arise from 
prior star formation activity in the halo (Ostriker & Thuan 1975) or the bulge. It was argued (Prantzos 2007), 
however, that while the Galactic halo reached a maximum metallicity of Z ~ 0.002 (^ 0.1 Z^), its mean 
metallicity is Z ~ 0.0006 (~ 0.03 Z.j|j) and its total mass (2x10® M^) is 20 times less than that of the disk 
(4.5x10^° M^). Koppen & Arimoto (1990) suggested a model in which the bulge evolves with a large yield 
(p « 0.034) and ends star formation by ejecting 1/10 of its mass as highly enriched gas with an abundance 
of 0.05 (~ 2.5 Z^) in a terminal wind which is captured by the proto-disk. Assuming that the proto-disk 
has a mass equal to that of the bulge, this mass transfer results in an initial disk abundance of 0.005 (^ 0.25 
Z.j;|j.) after mixing. Such a model is consistent with the metallicity distribution functions of both the bulge 
and disk (for more discussions, see P09). 

On the other hand, the p values derived from these models above are larger than the estimated value 
in the solar neighborhood: psoiN 0.7 Z^ = 0.014 (P09). We therefore need to establish models which not 
only match the metallicity distributions of the solar neighborhood but also give yields more consistent with 
observational data. 
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6.3. Inflow Models With Declining Rates 

In more realistic models, the disk does not evolve as a closed box and gas inflow or outflow (or both) is 
present. Inflow models have been attractive because of their capacity to reduce the number of low-metallicity 
stars born at early times, providing an elegant solution of the G dwarf problem, especially in view of the fact 
that gas accretion is expected to be a common phenomenon in the Universe (e.g., Prantzos 2007). There is 
indirect observational evidence for gas infall into the Galactic disk. The star-formation rate in the disk is 
about a few solar masses per year, while the total gas mass therein is ~ 5 x 10® M^. Thus, if there were no 
accretion of new material, our disk would run out of gas in a fraction of the Hubble time (Mo et al. 2010). 

The inflowing gas could originate from the gaseous halo, the bulge, the outer parts of the disk, or the 
accretion of small satellite galaxies. The disk itself may be considered as having formed completely from 
such inflowing gas or as having an initial mass (usually small compared to its final mass) with or without 
metals. Still more physically realistic models assume gas inflow with variable rates, and those with declining 
rates appear especially to provide in better agreement with recent observations (P09). In this section, two 
inflow models will be considered and the resulting metallicity distributions will be compared with that from 
this study. 


6.3.1. The Exponential Inflow Model (EIM) 

Some studies have shown that an exponentially decreasing infall rate with a long characteristic time 
scale ~ 7 Gyr can provide a reasonably good fit to the data. For example, GB04 demonstrated that such a 
model seems to reproduced the observed metallicity distributions of G and K dwarfs in the disk. 

For simplicity, some inflow models are based on a specific star formation rate given by 

= w;MG(t) (6) 

where t is time, MG(t) is the gas mass as a function of time and w is a constant. By defining a time-like 
variable u by 

u = wi, (7) 

Sommer-Larsen (1991) constructed a dynamical model with the assumption 

F(u) = Aw exp(—u) (8) 

in which F is the accretion rate of material from outside the system and A is a constant. If the initial mass 
and metallicity of the system are assumed to be zero and if the inflow gas has no metal content (Zp = 0) 
then, as u becomes very large, the expression for Z and dMs/dlog(Z) can simply be written as (P09): 

Z(u) = (^) (9) 


dMs 

dlog(Z) 


Z 2 

oc ( —) exp(—2Z/p) 
P 


( 10 ) 


The expression on the right side of the proportionality (10) is essentially the square of that corresponding 
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to the SCBM. The normalized distribution from the EIM for p = 0.025, which peaks at [Fe/H] = +0.1, is 
shown in Figure 11 (the red curve); it presents the metallicity distributions from this work (the histogram) 
and from the SCBM (the black curve) as well. Evidently, the EIM predicts fewer metal-poor M dwarfs and a 
steeper slope in the low metallicity tail than the SCBM, and is therefore in better agreement with the data. 
There is, however, an inconsistency between the value of p from this model and PsoiN, which is somewhat 
problematic. It should be remarked here that if our distribution had a peak at a lower metallicity, say, [Fe/H] 
« —0.1 dex, the model would result in a more reasonable yield. Accordingly, the accuracy of the metallicity 
distribution is critical for a comparison of GCE models. 



Fig. 11.— The comparison between the metallicity distributions of M dwarfs from this study (the histogram), 
the EIM (the red curve) and the SCBM (the black curve). 


6.3.2. Clayton’s Models 

Clayton (1985; 1987; 1988) has derived a series of models in which the inflow rate is parameterized by 

F(t) = ^MG(t) (11) 

t + to 

where k is an integer and to (or uo = icto) is arbitrary. By considering Equations (6) and (7), and assuming 
Zo = Zp = 0, expressions for Z and dMs/dlog(Z) in terms of u can be derived (Equations (8.38) and (8.39) 
in P09). The red curve in Figure 12 shows a distribution of this kind with k = 7 and p = 0.017, having 
a peak at [Fe/H] = +0.1 dex. The P value in this case is in better agreement with PsoiN- Although the 
model fits the low-metallicity tail of our distribution (the histogram) well compared with the SCBM (the 
black curve), it underestimates the number of supersolar-metallicity stars, and a modification of this model 
for the high-metallicity tail is needed. 

There have been other models such as merger models (Nagashima & Okamoto 2006) which resolve the 
G and K dwarf problems and can be tested using M-dwarf metallicity distributions in future. However, 
larger, more accurate data sets are required in order to discriminate among these models and find the more 
realistic one. 

The metallicity distribution of long-lived stars (i.e., G, K and M dwarfs) as well as the age-metallicity 
relation traced by the metallicity of these dwarfs are the most important observational constraints for GCE. 
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Fig. 12.— The comparison between the metallicity distributions of M dwarfs from this study (the black 
histogram), the Clayton’s model (the red curve) and the SCBM (the black curve). 


These low-mass stars belong to stellar populations of different ages and their metallicity distribution provides 
a fairly complete record of the evolutionary history of the Galaxy. However, CB04 pointed out that G dwarfs 
are sufficiently massive that they have begun to evolve away from the main sequence, and consequently, K 
dwarfs would make a cleaner sample for the local metal abundance distribution. For the same reason, we can 
expect that M dwarfs would provide an even better sample to represent the local metallicity distribution. 
M dwarfs are numerous, making up around half the total stellar mass of the Galaxy (RG97), and their 
metallicity distribution therefore offers a robust tool for studying the chemical evolution of the Milky Way. 
However, more precise metallicity calibrations are essential. 


7. SUMMARY 

Using a sample of 67 carefully selected dwarfs, we developed an optical-NIR photometric method to 
determine the metallicity of dwarf stars with spectral types between K6 to M6.5 in the metallicity range 
—0.73 to -1-0.3 dex. The calibration sample has two parts; the first part includes 18 M dwarfs in common 
proper pairs with an FGK star or early-type M dwarf of known metallicity and the second part contains 49 
dwarfs with metallicities obtained through an analysis of moderate-resolution spectra. Although the method 
may not be as accurate as moderate-to-high resolution spectroscopic techniques, they can be consistently 
applied to large numbers of stars, facilitating statistical investigations on metallicity distributions efficiently. 

We selected a large sample of around 1.3 million M dwarfs from the merged SDSS and 2MASS catalogs 
and corrected their gJK photometry for Galactic extinction based on a widely used dust map of the northern 
Galactic hemisphere. These stars meet the requirements for clean photometry and fall in color ranges 
necessary for removing possible giants as well as those ranges required for our metallicity calibration. We 
also applied a color cut needed for a photometric parallax to estimate stellar distances. Using the Z - height 
distribution, we found that the majority of M dwarfs in this large sample are found in the Galactic thin 
disk. By applying our calibration to these M dwarfs, the metallicity distribution of the local thin disk was 
determined and investigated. The metallicity distributions of subsamples with different Z - heights recovers 
the well known observation that there is a decrease in mean metallicity with Z - height, confirming yet again 
that older stars have on average lower metallicity and are farther from the Galactic plane. 

We examined the SGBM of GGE by comparing the distribution expected from this model with our 
M-dwarf distribution. A significant discrepancy between the model and our results in the low-metallicity tail 
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leads to the “M dwarf problem,” similar to the previously known G and K dwarf problems. More realistic 
models advanced to resolve the G and K dwarf problems may also be solutions to the M dwarf problem as 
well. We explored a few of these models such as the SCBM with pre-enrichment and two kinds of infall gas 
models with declining rates, and showed that these could more or less mitigate the M dwarf problem. 

In order to improve significantly these results and to determine a more reliable M-dwarf metallicity 
distribution, a more accurate calibration over a broader range of metallicities is required. The photometry 
of many M dwarfs of known metallicity is unusable because the images are saturated in SDSS images which 
accounts for their exclusion from our calibration sample and why only the g band was employed in the optical 
region. The acquisition of high-quality photometric griz data for a sample of M dwarfs with a variety of 
well-determined metallicities is already in progress. 
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Table 1. Observational Properties of 18 M dwarfs with metallicities determined by high-resolution spectroscopy 
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Table 2. Measured Properties of 18 M dwarfs with metallicities determined by high-resolution spectroscopy 
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et al. (2013a) in their calibration sample, based on the software package SME as mentioned above and a set of tuned lines from 
the SPOCS catalog (Valenti & Fischer 2005) through high-resolution observed spectra; WW05=Woolf and Wallerstein (2005) 
through high-resolution observed spectra 



Table 3. Astrometry® and Extinction-Corrected Photometry of 49 dwarf stars in the calibration sample 
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KIC 5252367 283.1939 40.4956 16.192 0.004 un.sat 11.654 0.027 11.022 0.030 10.847 0.018 
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Table 4. Spectral Type and Metallicity of 49 dwarf stars in the calibration sample 


Name 

Spectral Type 

Spectral Type Ref. 

[Fe/H] 

[Fe/H]Ref. 

LSPM J0856+1239 

M6 

N14 

+0.07 

± 

0.12 

N14 

LSPM J1005+1703 

M6 

N14 

+0.06 

± 

0.12 

N14 

LSPM J1530+0926 

M6 

N14 

+0.09 

± 

0.12 

N14 

LSPM J2012+0112 

M6 

N14 

+0.06 

± 

0.12 

N14 

LSPM J0738+1829 

M6 

N14 

+0.30 

± 

0.12 

N14 

LSPM J1002+4827 

M6 

N14 

+0.28 

± 

0.12 

N14 

LSPM J1432+0811 

M6 

N14 

+0.24 

± 

0.12 

N14 

LSPM J1631+4051 

M6 

N14 

+0.30 

± 

0.12 

N14 

LSPM J0738+4925 

M6 

N14 

+0.17 

± 

0.13 

N14 

LSPM J1239+0410 

M6 

N14 

+0.30 

± 

0.12 

N14 

LSPM J1309+2859 

M5 

N14 

+0.08 

± 

0.13 

N14 

LSPM J1021+0804 

M5 

N14 

+0.10 

± 

0.12 

N14 

LSPM J1352+6649 

M5 

N14 

+0.15 

± 

0.12 

N14 

LSPM J0918+6037W 

M5 

N14 

+0.10 

± 

0.12 

N14 

LSPM JOOOl+0659 

M5 

N14 

+0.10 

± 

0.12 

N14 

LSPM J0237+0021 

M5 

N14 

+0.03 

± 

0.12 

N14 

LSPM J0917+5825 

M5 

N14 

+0.08 

± 

0.12 

N14 

LSPM J1031+5705 

M5 

N14 

+0.01 

± 

0.13 

N14 

LSPM J1112+0338 

M5 

N14 

+0.04 

± 

0.13 

N14 

LSPM J1419+0254 

M5 

N14 

+0.04 

± 

0.13 

N14 

LSPM J0024+2626 

M5 

N14 

+0.29 

± 

0.12 

N14 

LSPM J1348+0429 

M5 

N14 

+0.27 

± 

0.12 

N14 

LSPM JllOl+0300 

M5 

N14 

-0.30 

± 

0.13 

N14 

NLTT 8870 

M5 

M14 

+0.11 

± 

0.10 

M14-13a Empir.Calib 

LSPM J0959+4712 

M4 

N14 

+0.01 

± 

0.12 

N14 

LSPM J1148+5305 

M4 

N14 

-0.01 

± 

0.13 

N14 

LSPM JlllO+4757 

M4 

N14 

-0.09 

± 

0.13 

N14 

LSPM J0958+0558 

M4 

N14 

+0.03 

± 

0.12 

N14 

LSPM J1316+2752 

M4 

N14 

-0.20 

± 

0.14 

N14 

LSPM J1709+3909 

M4 

N14 

-0.23 

± 

0.14 

N14 

LSPM J1240+1946 

M4 

N14 

-0.14 

± 

0.14 

N14 

LSPM J0920+0322 

M4 

N14 

-0.28 

± 

0.14 

N14 

LSPM J1345+2852 

M3.5 

D13/R10-13 

-0.13 

± 

0.12 

D13 

LP 655-23 

M3 

M14 

+0.03 

± 

0.10 

M14 

LSPM J1314+4846 

M3 

N14 

-0.13 

± 

0.12 

N14 

LSPM J1038+4831 

M3 

N14 

+0.24 

± 

0.12 

N14 

LSPM J2148+0126 

M3 

N14 

+028 

± 

0.12 

N14 

KIC 5252367 

M2 

P12/R10-13 

-0.18 

± 

0.13 

M13b/Non-KOI sample 

KIC 6183511 

Ml.5 

P12/R10-13 

-0.15 

± 

0.06 

M13b/KOI sample 

KIC 3426367 

Ml.5 

P12/R10-13 

-0.15 

± 

0.13 

M13b/Non-KOI sample 

KIC 3935942 

Ml 

P12/R10-13 

+0.04 

± 

0.05 

M13b/Non-KOI sample 

KIC 4543236 

Ml 

P12/R10-13 

-0.19 

± 

0.06 

M13b/Non-KOI sample 

KIC 4243354 

MO 

P12/R10-13 

-0.33 

± 

0.12 

M13b/Non-KOI sample 

KIC 5513769 

MO 

P12/R10-13 

-0.05 

± 

0.09 

M13b/Non-KOI sample 

KIC 4725681 

MO 

P12/R10-13 

-0.01 

± 

0.07 

M13b/KOI sample 

KIC 4139816 

K7.5 

P12 

-0.51 

± 

0.07 

M13b/KOI sample 

2MASSJ20515725“ 

K7 

SSPP/ST 

-0.58 

± 

0.05 

SSPP/M 

KIC 4543619 

K6.5 

MQ14 

-0.55 

± 

0.12 

M13b/Non-KOI sample 

KIC 2850521 

K6 

P12 

-0.68 

± 

0.11 

M13b/Non-KOI sample 


“2MASSJ20515725-0111317 


Note. — Source of Spectral types: M14=Mann et al. (2014), N14=Newton et al. (2014), 
D13/R10—13=Temperature taken from Deshpande et al. (2013) and spectral type determined 
by the temperature-spectral relations of RIO and R13, P12/R10—13=Temperatures taken from 
Pinsonneault et al. (2012) and spectral type determined by the temperature-spectral relation 
of RIO and R13, P12= Spectral type estimated by temperatures taken from Effective tem¬ 
perature scale for KIC stars (Pinsonneault et al. 2012), SSPP/ST= Spectral type estimated 
by temperature taken from the Sloan Extension for Galactic Understanding and Exploration 
(SEGUE) Stellar Parameter Pipeline (SSPP, Lee et al. 2008), and MQ14=Spectral type esti- 
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mated by temperature taken from Rotation periods of Kepler main sequence stars (McQuillan 
2014) Source of [Fe/H]: N14=Based on an empirical metallicity calibration of moderate- 
resolution, NIR spectra in Newton et al. (2014), M13b/KOI sample=Tlie sample taken from 
Table 1 in Mann et al. (2013b) with metallicities based on a modified empirical calibration 
of moderate-resolution, optical spectra in Mann et al. (2013a;2013b), M13b/Non-KOI sam- 
ple= The sample taken from Table 2 in Mann et al. (2013b) with metallicities based on the 
weighted means of J-, H-, and 77-band calibrations described in Mann et al. (2013a), and 
SSPP/M= Based on multiple approaches of SSPP; M14-13a Empir.Calib=Mann et al. (2014), 
based on the empirically spectroscopic calibration through moderate-resolution observed spec¬ 
tra; Rob07=Robinson et al. (2007) through moderate-resolution observed spectra 
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